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I. INTRODUCTION 

Spin properties of few electron quantum dots have 
been recently extensively studied, in hope that a local- 
ized spin can serve as a qubit, a central building block 
of a quantum computer. Spin, compared to orbital de- 
grees of freedom, was anticipated to have much longer 
coherence time. Fast experimental progress during last 
few years supported this assumption - long electron spin 
relaxation^iii^ and dephasing times^ii have been mea- 
sured in quantum dots. 

If a quantum dot electron spin is to realize a qubit, 
DiVincenzo's criteria have to be fulfilled:^ (i) The exis- 
tence of a qubit — the two states of spin naturally encode 
the information bit. (ii) The initialization of the qubit is 
also straightforward - at a finite static magnetic field and 
small temperature it is enough to let the system relax into 
the ground state spontaneously, (iii) Due to the isolation 
of the spin from the environment, the qubit readout is not 
that easy, but can be now considered experimentally mas- 
tered, using spin-to-charge conversion schemesi^iiS (iv) 
Concerning coherent manipulation, a very important step 
forward is a recent demonstration of magnetically driven 
Rabi oscillations ?ii Thus, all basic ingredients have been 
shown to work at the proof of principle level, and the ef- 
fort now is aimed at their integration, with the final goal 
of a (v) scalable qubit design. Connected with the last 
two points, namely, if dots in an array can be addressed 
individually, the manipulating fields have to be produced 
locally, nearby the particular dot being manipulated (so 
called fields generated on-chip). 

If the spin is manipulated by a magnetic field which 
is produced locally (by an oscillating current in a wire 
nearby the dot) , it is inevitably accompanied by an oscil- 
lating electric field. This electric field is due to an imper- 
fect screening of the dot from the circuitry; the electric 
field due to a changing magnetic field, V x E = — 9tB 
is negligible, The electric field strongly disturbs the or- 
bital part of the electron wavefunction and, if spin-orbit 
coupling is present, also couples to the spin - one way 



or the other, it limits the strength of the magnetic field 
usable for the manipulation (in Ref. TTthis limit was 1.9 
mT) and thus limits speed of the operation (the maximal 
achievable Rabi frequency). 

On the other hand, the electric field does not have to 
be viewed as an enemy - through the spin-orbit interac- 
tion it induces the very same spin oscillation o^^d^ as the 
magnetic field. From practical point of view, the elec- 
tric field is preferred, since it is much easier to control 
than the magnetic field. The possibility of electrically 
induced spin oscillations is eagerly pursued experimen- 
tally. The disadvantage of the electric field is that the 
Rabi frequency is expected to depend on dot parameters, 
since the coupling is due to material dependent spin-orbit 
interactions. (In the case of the magnetic field, the fre- 
quency is given only by the field strength.) 

It is thus an important issue to compare the effective- 
ness of oscillating electric and magnetic fields in inducing 
Rabi spin oscillations. Namely, how large fields are re- 
quired to induce Rabi oscillation with certain frequency 
and how the frequency depends on the parameters of the 
dot (consequently, how stable it is against fluctuations of 
these parameters). This is where this paper aims - we 
quantify the dipolar electric and magnetic couplings in 
spin resonance of a single electron confined in a quantum 
dot. By spin resonance we mean that a static magnetic 
field is applied, whereby the ground orbital state is split 
by the Zeeman energy. Oscillations between the two split 
states are induced by oscillating magnetic and electric 
fields if the field frequency equals to the Zeeman energy. 

It was already showed theoretically in single dots that 
due to the presence of spin-orbit interactions the elec- 
tric field is indeed effective in inducing spin resonance 
For a typical lateral single GaAs quantum dot, in static 
magnetic field of 1 Tesla, an oscillating electric field 
of 10"^ V/m is as effective as the oscillating magnetic 
field of I mT. We widen the analysis on the experimen- 
tally relevant double dot setup. Our main result is that 
the anisotropy of the spin-orbit interactions allows for a 
strong control over the electrical field efficiency in spin 
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manipulations. Our findings provide guide for dot con- 
figurations for two possible strategies: If a local electric 
field is chosen for the spin manipulation, we show how 
its efficiency can be maximized. On the other hand, if 
a magnetic field is chosen, the coupling due to the ac- 
companying electric field is unwanted - we show that 
it can be suppressed by: (i) lowering the magnitude of 
the static magnetic field, (ii) special orientation of the 
static magnetic field. Particularly, in an easy passage 
configuration)^ the otherwise most effective electric field 
component is completely blocked in disturbing the spin. 

In addition to spin resonance, we also study the electri- 
cally and magnetically induced orbital resonance, where 
the resonant states are the two lowest orbital states with 
the same spin. A qubit represented by these two states 
is called a charge qubit. The study is motivated by an 
observation, that in the presence of spin-orbit interac- 
tions, an analogue to electrically induced spin resonance 
should exist. Namely, the magnetic field should induce 
oscillation between spin alike states. We show that this 
is indeed true, however, for realistic values the magnetic 
field is much less effective compared to the electric field. 

We use realistic parameters for electrically defined sin- 
gle and coupled dots in [001] grown GaAs heterostruc- 
ture. We treat the problem numerically by exact diago- 
nalization of the full electron Hamiltonian. However, for 
all results we provide explanatory analytical arguments 
based on an effective spin-orbit Hamiltonian and the de- 
generate perturbation theory. Our model incorporates 
the electron relaxation and decoherence rates caused by 
acoustic phonons in a realistic way; the rates we use have 
been found to be in a very good agreement with the ex- 
perimental data for magnetic fields above 1 Tesla both in 
singl^iii. and double dotsj^ 

The paper is organized as follows: In Sec. II we 
describe the model of the electron in the dissipative 
phonon environment under oscillating electric and mag- 
netic fields. In Sec. Ill we derive an effective spin-orbit 
Hamiltonian which allows symmetry analysis of the prob- 
lem. With this Hamiltonian we evaluate the matrix el- 
ements of oscillating magnetic and electric fields for the 
case of spin (Sec. IV) and orbital (Sec. V) resonance. In 
Sec. VI we describe the system in the steady state, where 
we show how to obtain the Rabi frequency and decoher- 
ence from a steady state measurement. 

II. MODEL 

Consider a single electron in a double quantum dot 
formed in a two dimensional electron gas in a (001) plane 
of a GaAs / AlGaAs heterostructure. The effective Hamil- 
tonian is 

H = Ho + Hbr + Hd + Hds, (1) 

where 

Ho=T + Vc + Hz. (2) 
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FIG. 1: The orientation of the potential dot minima (denoted 
as the two circles) with respect to the crystallographic axes 
{x = [100] and y = [010]) is defined by the angle S. The 
orientation of the inplane magnetic field is given by the angle 
7- 

The kinetic energy is T = fik^/2m with the effective elec- 
tron mass m and kinetic momentum Kk — —iK^ . The 
double quantum dot is described here by the confinement 

Vc{v) = imwo'min{(r - d)^, (r + d)^}, (3) 

representing two alike potential minima of a parabolic 
shape, centered at ±d. The in-plane orientation of the 
double dot with respect to the crystallographic axes x 
and y is defined by 5, the angle between d and x. A 
single dot with the confinement energy = hioo, and 
the confinement length Iq ~ y/h/mujQ, is defined by the 
limit d ~ 0. Alternatively to giving the interdot distance 
d, the double dot can be characterized by tunneling en- 
ergy SEt equal to half of the difference of the energies 
of the two lowest orbital statesJ^ The electron feels an 
in-plane magnetic field B whose orbital effects can be ne- 
glected, for fields lower than ~ 10 T. The Zeeman term 
is Hz — /icr.B, where fi = {g/2)fj,B is the renormalized 
magneton, g is the conduction band g- factor, /x^ is the 
Bohr magneton, and cr are the Pauli matrices. The spin 
quantization axis is defined by the direction of the mag- 
netic field. The angle between B and x is denoted as 7. 
The geometry is summarized in Fig. [TJ 

The spin-orbit coupling in our confined system 
is described by three terms.— The Bychkov-Rashba 
Hamiltonianjiiii^ 

Hbr = — -, — (o-xky - Uyk^) , (4) 
2mlBR 

is present due to the heterostructure asymmetry, while 
the linear and cubic Dresselhaus HamiltonianS )^^'^" 

Hd = T, — r- {-o-xkx + cTyky) , (5) 
ZmiD 

Hd3 = 7c {<Jxkxkl - CTyky kl) , (6) 
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are due to the lack of the bulk inversion symmetry. 

In our numerical calculations we use bulk GaAs mate- 
rial parameters: m — 0.067 me, g = —0.44, and 7c — 27.5 
eVA^. For the coupling of the linear spin-orbit terms we 
choose Ibr = 1-8 /im, and lu = 1.3 fiia, the values used 
to fit a recent experiment For the confinement length 
we take Iq = 30 nm, corresponding to the confinement 
energy Eq — 1.2 meV. 

We now describe the influence of the phonon environ- 
ment as well as of the oscillating electric and magnetic 
fields. The phonon environment leads to the relaxation 
and decoherence expressed, in the Markov and Born ap- 
proximations, by the time derivative of the diagonal and 
off-diagonal elements of the reduced density matrix of the 
electron, (upperscript "ph" stands for phonons, to 
discriminate from other contributions to the time deriva- 
tive which appear later) 

df^ Pa = -^"^^ikPu +^2rkiPkk, (7a) 

k k 

df^Pij = -^i^ik+TJk)p^J = -lijPtj- (7b) 

k 

Here 2rij is the relaxation rate from the electron state 
i to j due to the piezoelectric and deformation potential 
interactions of the electron with acoustic phonons. There 
is no additional phonon channel for the decoherence jij 
apart from the relaxation, since the phonon density of 
states vanishes for zero phonon energy, Tu — 0. We 
do not consider non-phonon mechanisms of dephasing, 
which are important at low (sub Tesla) magnetic fields. 
To allow for a finite temperature one can suppose a de- 
tailed balance: Tji — rFy , where r = eyi^{—hwij/kBT). 
In the calculations below, we consider temperature much 
lower than the orbital excitation energy. For example, 
the experiment Ref . fill was done at temperature 100 mK, 
corresponding to ~ 0.01 meV, while a typical excitation 
energy of the used quantum dot was about 1 meV. In 
this limit a transition into a higher orbital level has a 
negligible rate. 

In addition to phonons, the electron is subject to os- 
cillating electric and magnetic fields, which contribute 
through the following Hamiltonian: 

= [eC.r + /iS.cr] cosujt = Mlcosujt. (8) 

Only the in-plane components of the oscillating electric 
field are relevant. The oscillating magnetic field is per- 
pendicular to the plane, B cx z, simulating the condi- 
tions in the experiment;^ In the numerical calculations 
we set E ~ 1000 V/m as a realistic guess for the ex- 
perimental setup^^ and B — \ mT, a typical value from 
the experiment ii We suppose that frequency lo is close 
to the energy difference of a given pair of states - res- 
onant states - denoted by indexes a and &, such that 
Lj w LJiia = {Eb — Ea)/h > 0. In the rotating wave 
approximatiouf^i that we adopt, the oscillating field in- 
fluences only the two resonant states, contributing to the 



time derivative of the density matrix: (upperscript "of" 
stands for the oscillating field) 

dfpaa = -dfpbb = ^flbaPabe'''' " ]^^abPbae-''^\^a.) 
dfpab = -\{pbb- Paa)^abe~''^\ (9b) 

where A = ujba — w is the detuning from the resonance. 

The time evolution of the density matrix, given by 
Eqs. (O and © can be easily solved if one neglects all 
other states but the two resonant;^ Such approxima- 
tion makes sense if the electron can not escape from the 
two state subspace. Roughly speaking, the effective rate 
out of the subspace must be much smaller than rates for 
transitions restoring the electron in. This, for example, 
means that the ground state must be one of the reso- 
nant states, which is the case here. Another interesting 
counterexample is optical shelving,-24 whereby the elec- 
tron can be trapped in an intermediate dark state. There 
are parameter values for our model where the three lowest 
electron states can realize such a scheme, but we do not 
discuss it in this article. We work in the regime where 
the two level approximation is justified, as follows also 
from our numerical results. The validity of the two level 
approximation also implies that the decoherence rate is 
given by the relaxation only, 

lab = Iba = Tfca + Tafc = Tbai^ + t), (10) 

a fact that we will use later. 

Suppose now the electron is in the ground state ini- 
tially. After the resonant field is turned on, the popula- 
tions of the two resonant states start to oscillate, mean- 
ing, after a certain time the electron will be in the ex- 
cited state, then comes back to the ground state and so 
on. Since these Rabi oscillations are coherent, they can 
realize a single qubit rotation, one of the basic building 
blocks of a quantum computation. The time after which 
the populations switch is proportional to the inverse of 
the frequency of the Rabi oscillations (Rabi frequency) fl. 
A larger Rabi frequency then means a faster single qubit 
operation. To better assess the suitability for quantum 
computation, one has to take into account also the decay 
of the Rabi oscillations which is due to the decoherence. 
In our model, the magnitude of the oscillations decays 
exponentially with the rate roughly proportional to the 
decoherence rate jba- Therefore, to minimize the error 
in a single qubit operation, it is desirable to maximize 
the ratio fl/^bai which quantifies how many single qubit 
operations one can do during the decoherence time. We 
note that from the observed decaying Rabi oscillations in 
the time domainj^iii both il and ^ba can be extracted. 

As the last here we note that in the two resonant states 
approximation there are three important rates, decoher- 
ence jba, detuning A and the field matrix element \^lab\- 
If the the last one is not dominant, then either jab is 
large and the damping is too strong to observe Rabi os- 
cillations or A is large and the magnitude of the oscilla- 
tions is small^ - both cases are not of interest here. We 
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FIG. 2: (Color online) The lowest part of the energy spec- 
trum of Hamiltonian Hq, Eq. ([2]), in zero magnetic field as a 
function of the tunneling energy (SEt) in the units of the con- 
finement energy Eq . Each eigenfunction belongs to one of the 
four symmetry classes of C2v, which are denoted by different 
color/type of the line. Spin indexes are omitted. 
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Higher order terms and a constant were omitted in Hi. 
The curly brackets denote the anticommutator, while hi 
is an effective spin-orbit vector specified below. 

For the following discussion the symmetries of the 
terms in Eq. p4p are important. First, each term has 

(2) 

a definite time reversal symmetry: is antisymmet- 

ric, while the other terms are time reversal symmetric. 
Second, the spatial symmetry of a particular term is de- 
fined by a combination of x, y, k^, and ky it contains. To 
exploit the spatial symmetry of the confinement, Eq. ([3]), 
we rotate the (originally crystallographic axes) coordi- 
nates such that the new x lies along d. The coordinates 
change according to 



consider the case when the field matrix element is indeed 
dominant. It holds then that the matrix element equals 
the Rabi frequency, 

n = \nba\, (11) 

and is therefore of crucial importance. In the next we 
analyze in detail how the field matrix element due to 
electric and magnetic oscillating fields depends on system 
parameters. To simplify the analysis of the spin-orbit 
influence, we begin with a derivation of an effective spin- 
orbit Hamiltonian. 



III. EFFECTIVE SPIN-ORBIT HAMILTONIAN 

It is useful to remove the linear spin-orbit terms in 
Eq. ([T]) by applying a unitary transformation^^ leading 
to a new Hamiltonian 

H' ^ e^Her^ ^ Hq + Hi, (12) 

where 

S=7^{y(^x-xay)- ^{xa^ -yay) (13) 

is a transformation matrix and 

Hi = Hd3 + hI^,! + Hf + H^^l (14) 

is an effective spin-orbit interaction. In addition to the 
cubic Dresselhaus term Hj^^ , Hi comprises the following 



a; — > a; cos (5 — y sin (5, — > y cos (5 + a: sin (5, (18) 

and similarly for kx and ky . The rotation leaves Eq. (jlSp 
unchanged. In Eq. (|16p the effective linear spin-orbit 
couplings and h\ acquire the following form: 

hi = ls],cosi-f~d)-la'sm{-f + S), (19) 
hi = lg]^sm{j-6)-l]^^cos{-f + 6). (20) 

It is important that these couplings can be selectively 
turned to zero by orienting the static magnetic field B 
in a certain direction (7) dependent on the orientation of 
the double dot (S) . The result of the rotation in Eq. pT]) 
is not presented here; we will discuss only its relevant 
parts. 

We can obtain analytical results in reasonable quanti- 
tative agreement with the numerics in the lowest order 
degenerate perturbation theory by exploiting the symme- 
tries of the problem. The orbital eigenfunctions of Hq, 
Eq. ([2]), in an in-plane magnetic field form a representa- 
tion of symmetry group. There are four possible 
symmetry classes which transform upon inversions along 
(rotated axes) x and y as l,x,xy, and y, respectively. 
A relevant part of the double dot spectrum is shown in 
Fig. [21 Several eigenstates are labeled by F with indexes, 
where the bottom index denotes the spatial symmetry of 
the state (four symmetry groups), while upper indexes 
labels states within the symmetry group - this notation 
was introduced in Ref. llSl . In further the two lowest or- 
bital states will play the most important role: ground 
state F™ is symmetric both in x and y (often denoted 
as the bonding molecular orbital), and first excited or- 
bital state F2° is antisymmetric in x and symmetric in y 
(antibonding). 
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If a magnetic field is applied, each line in Fig. [2] splits 
into two by the Zeeman term lifting the degeneracy. As- 
suming negative g-factor and positive i?, a spin down 
state (denoted by |) has higher energy than a spin up 
state (t). A further important consequence of a finite 
Zeeman energy is the anti-crossing of states and Fj^ , 
influence of which we take into account using the degener- 
ate perturbation theory. Exact eigenfunctions (denoted 
by an overline) of the Hamiltonian H' can be written as 
a combination of the solutions of Hq (denoted by F as in 
Fig. [2]): the three lowest states, in the lowest order of the 
degenerate perturbation theory, are 



-00 



=00 



-10 

2T 



T^OO 



n 



q;F2| 



e;10 



■ 2i 



Efl El] 



4T 



10\ 



/3*F?? 



tpOO 
^2T 



-a 



(21) 
(22) 
(23) 



The dots denote the rest of an infinite sum through the 
eigenfunctions of Hq. The anti-crossing is described by 
coefRcients 

13 = Arg(c(5i;)sin[arctan(|4c/(5£;|)/2], (24) 
a - yT~W, (25) 
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which depend on the energy difference 5E 
and the coupling c = (Fji-ffiF™) between the unper- 
turbed crossing states. 

From the above expression for coupling c it follows that 
the anti-crossing is caused by the part of Hi with the spa- 
tial symmetry of x, which is the symmetry of F2|. After 
the rotation, Eq. (jlSp . the only term with x symmetry 

in Eqs. ((I5l)-(IT71) is the first term in H^^\ Therefore, 
by orienting the magnetic field such that = 0, one 
can turn the anti-crossing into a crossing, /3 = 0. Note 
that also -ffca contains a term of x symmetry - this does 
not hinder to achieve = 0, but only slightly shifts the 
required position of the magnetic field. Changing the 
anti-crossing into a crossing has profound conse que nces 
on the spin relaxation time, as was found in Ref. As 
we will see here, it is similarly important also for the 
electrically induced spin resonance. 

We finish this section with a note about other formula- 
tions of the unitary transformation Eq. . It was first 
used in the context of quantum dots in Ref 25. There 
they neglected the cubic Dresselhaus term, but kept the 
correction of the third order in the spin-orbit couplings, 
which in our notation would be 



(2)i 



(26) 



This term, which we neglected, together with Eg. p^ . 
were there interpreted as a vector potential of a spin- 
orbit originated magnetic field. 



If the potential is harmonic (d = in out model) , the 
unitary transformation can be generalized to remove ex- 
plicitly also the lowest order mixed Zeeman-spin-orbit 

(2) 

term at the expense of parameters of the Hamilto- 
nian (like mass) becoming spin dependent. However, this 
possibility is specific to the potential form and nothing 
can be done with the cubic Dresselhaus term. 

An elegant form of the unitary transformation together 
with the perturbation theory is worked out in Ref. [T^ . 
where an effective Hamiltonian for a set of degenerate 
states is derived in a compact form using an inverse of Li- 
ouville superoperator. However the inverse is not known 
for any other potential than harmonic if the Zeeman term 
is present and is not known at all for the cubic Dressel- 
haus term. 

The effective Hamiltonian presented here is indepen- 
dent of the confinement potential form and reveals the 
symmetry of the spin-dependent perturbations. In a 
symmetric potential, such as our double dot, simply by 
inspecting the symmetry of the terms allows to identify 
the term responsible for certain process (such as spin re- 
laxation, or electrically induced transition). Formulas in 
Eqs. (fT5|) - P?)l hold also if an out of plane component of 
the magnetic field Bz is present, provided that (i) the op- 
erator k is promoted to contain also the vector potential 
of this component k — ifiV 4- eA(_B^) and (ii) there is 

(2) 

an additional contribution to , proportional to B^, 
see Ref. 26 for its form. 

Using the effective spin-orbit Hamiltonian and the ap- 
proximations for the eigenstates, we now quantify indi- 
vidual contributions of oscillating fields to the matrix el- 
ement flba ■ We will show where these contributions orig- 
inate and how they can be used to control the electron 
spin and orbital resonance. 



IV. MATRIX ELEMENTS: SPIN RESONANCE 



As the spin resonance we denote a situation, when the 
two resonant states are the ground state T^^ and its Zee- 
man split counterpart . In this case we tag the matrix 
element by subscript "spin" , 



spin 



(27) 



where the upperscript F will stand for the particular part 
of f2 we consider. But before dealing with the specific 
oscillating field Hamiltonian, we remind the Van Vleck 
cancellation that occurs in Eq. (P7)l due to the time re- 
versal symmetry. Consider a general Hermitian operator 
O. Let us write the matrix element in the following form: 



(r^loiF-) 



(f;°|6|f°;) + <5o. 



(28) 



where 60 is due to spin-orbit corrections. If the first 
term is nonzero, that is the unperturbed states are cou- 
pled by O, SO can be usually neglected. If the first term 
vanishes, and we are away from the anti-crossing, /? ^ 1, 
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the time inversion symmetry gives an important infor- 
mation about the matrix element 60. Indeed, if O has a 
definite time reversal symmetry, T{0) = f (— 1) when be- 
ing symmetric (antisymmetric), using Eqs. (PT|) - (P^ for 
the matrix element in the lowest order in Hi we get^^'^* 



I T{Hi)T{d) \ 



(29) 



where i denotes the symmetry class, j denotes, for 
brevity, both upper orbital indexes, and a denotes the 
spin. In this lowest order, the contributions from the 
constituents of Hi are additive and can be considered 
separately. Therefore the first order contributions of the 
terms with the same time reversal symmetry as O [that is 
iiT{Hi)T{0) — 1] will be suppressed by a factor of order 
of ijlB/Eq, compared to matrix elements such as Eq. ([27|) . 
but between states with different spatial indexes. Near 
the anti-crossing the terms containing coefficients a and 
/3 dominate other terms in Eqs. and the matrix 

elements are then proportional to these coefhcients - the 
Van Vleck cancellation does not occur. 

These general results can be applied to the spin reso- 
nance due to magnetic and electric fields. The oscillat- 
ing magnetic field [Ml = fiBzCz] couples the unperturbed 
states: 



(30) 



so that we can neglect the spin-orbit contribution to the 
matrix element, dfl. 

On the other hand, the electric field dipole operator 
{hCl = e£.r) does not couple the unperturbed states. As 
Cl is now time reversal symmetric, the contributions of 

(2) 

all terms in Hi but H^ are suppressed. For the electric 
field along the rotated x axis the matrix element at the 
anti-crossing is 



(31) 



Away from the anti-crossing. 



spin 



2{Ei - El^) 



(Ei - Ef^Y 



(32) 



The spatial symmetry (here x) of the dipole operator se- 
lects only eigenfunctions of symmetry x in the sum. Only 
H^ , Eq. (|16p , contains a term of x symmetry, propor- 
tional to hf . In the above sum each state j (with energy 
E2) contributes proportionally to its dipole matrix ele- 
ment Xj. To get the analytical result close to numerics 
one needs to include the two lowest eigenfunctions in the 
sum in Eq. (|5^ . 

If the electric field is along the rotated y axis, the anti- 
crossing does not influence the matrix element, since y 
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FIG. 3: Calculated matrix elements between the resonant 
states due to magnetic and electric oscillating fields. The two 
upper panels, (a) and (b) , show the matrix elements f^spin for 
the spin resonance, while the two lower panels show orbital 
resonance elements f2orb- On the left, in (a) and (c) the ele- 
ments are functions of the static magnetic field, with a fixed 
tunneling energy of 20% of the confinement energy. On the 
right, in (b) and (d) the elements are functions of the tunnel- 
ing energy at a fixed magnetic field _B = 1 T. The dots are 
oriented along [100], while the static magnetic field lies along 
[010]. 



dipole operator of the electric field does not couple the 
ground and anti-crossing states. Then, an analogous ex- 
pression to Eq. p2p holds at (up to factor a multiplying 
some terms in the sum) or away from the anti-crossing: 



spin 



-e8yhy^,BJ2\Y,\ 



2{El - Ef') 



{Ei - Ef^Y - 4{fiBY 



(33) 



Here it is enough to include just the lowest eigenfunction 
of y symmetry in the sum. The dipole elements and en- 
ergy differences, computed by approximating the unper- 
turbed functions T by symmetrized single dot orbitalsf^^ 
are summarized in Tab. [H 

Fully numerical results for the matrix elements as a 
function of the magnetic field are shown in Fig. [3^. The 
matrix element of the magnetic field is constant, up to 
a narrow region of suppression due to a, since it de- 
pends only on the strength of the oscillating magnetic 
field, Eq. p0|) . The matrix elements of the electric field 
[Eqs. ([5^ and ([55)1 ] are proportional to the Zeeman en- 
ergy fiB - the spin resonance is more sensitive to elec- 
trical disturbance as the magnetic field grows, while at 
zero magnetic field the electric field is ineffective. At the 
anti-crossing, ^fp^^ is strongly enhanced (by two orders 

of magnitude) and described by Eq. (PT|) , while fifpi„ de- 
velops a small dip similar to ^^pin- 

It can be seen in Fig. where the matrix elements 
are functions of the tunneling energy, that the spin res- 
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TABLE I: Analytical approximations for the dipole matrix 
elements and energy differences. For each quantity the defi- 
nition, unit, expression, and limits for small and large inter- 
dot distances are given. In some cases the expression is too 
lengthy and only the asymptotics are given. The expression 
for 5Et is given in Ref. fl5|. The interdot distance measured 
in the units of the confinement length is used, T> = d/lo. 



onance is much more sensitive to the electric field along 
the double dots x axis than to a perpendicular field. This 
difference is strengthened at the anti-crossing. Only in 
the truly single dot case (d = or d = 00), the electric 
field influence is isotropic. We can also conclude from the 
single dot values that the matrix elements of magnetic 
field of 1 mT and electric field of 10'^ V/m are compa- 
rable in magnitude in the static magnetic field of order 
of Tesla. This means that in the experiment fii where no 
electrically induced signal was observed, the electric field 
is likely considerably lower than the estimated value of 
10^ V/m. 

Similarly to the spin relaxation rates^^i^^ the matrix 
element of the resonant electric field is highly anisotropic. 
The possible control over the resonance is demonstrated 
in Fig. 3^, where the matrix elements are shown as func- 
tions of the orientation of the static magnetic field. The 
magnetic field matrix element is independent on 7, as fol- 
lows from Eq. ([50]) . The electric field matrix elements are 
anisotropic, with the dependence given by the effective 
spin-orbit couplings and h\. By proper orientation of 
the static magnetic field it is thus possible to turn off the 
contribution due to the electric field pointed along a given 
direction. In particular, the electric field along x is not ef- 
fective {hi = 0) at 7 = arctan(/£)/ZB/i) « 38°. The elec- 
tric field along y is ineffective if 7 = arctan(ZB/j/?£)) 
58°, since here h\ — 0. These conditions were obtained 
from Eqs. HH) and ^ by putting 5 = (the dots 
oriented along [100]). Different orientation of the dots 
changes the conditions for the effective spin-orbit cou- 
plings to be zero. For example, in Fig. [JJd, the dots are 
oriented along [110], that is (5 = 45° and the effective 
coupHngs /if and h\ are zero at 7 = 45° and 135°, re- 
spectively, independent on the spin-orbit parameters. If 
the electric field points along a general direction, it is still 
possible to turn off the matrix element by properly ori- 
enting the magnetic field. However, in a general case the 
desired position of the magnetic field is defined not only 
by the effective couplings /if and but by all terms in 




y [deg] y [deg] 



FIG. 4: Calculated matrix elements for the spin [upper two 
panels (a) and (b)] and the orbital [lower two panels (c) and 
(d)] resonance due to oscillating magnetic and electric fields 
as functions of 7, the orientation of the static magnetic field, 
B = 1 T. The tunneling energy is 20% of the confinement 
energy. On the left, in (a) and (c) the dots are oriented along 
[100], that is 5 = 0. On the right, in (b) and (d) the dots are 
oriented along [110], 5 = 45°. 



Eqs. (1321)- iSl)- 

In the easy passage configuration, defined by = 0, 
the spin relaxation time does not suffer a drastic suppres- 
sion due to the anti-crossing, as was shown in Ref. [ij We 
see that in addition to that here the spin resonance is in- 
sensitive to otherwise most effective electric field compo- 
nent - along X. Such electric fields are inevitably present 
if the spin qubit is manipulated by an on-chip generated 
magnetic fieldJ^ On the other hand, on-chip manipula- 
tions seem inevitable in a scalable system, where it must 
be possible to address the qubits selectively. The easy 
passage configuration thus protects the spin against the 
electric field and provides a stable Rabi frequency over a 
wide range of parameters values, if the qubit is manipu- 
lated by an oscillating magnetic field. 



V. MATRIX ELEMENTS: ORBITAL 
RESONANCE. 

As the orbital resonance we consider the case when the 
resonant states are the two lowest orbital states and the 
matrix element is 

n^^,,, = {t^,mff,). (34) 

A similar suppression as in Eq. (|29p takes place also now, 
if the operator O acts only in the spin subspace (that is, 
it is the Zeeman term) . This suppression again favors the 

contribution due to iJ^ ' compared to the rest of Hi . If 
the anti-crossing dominates, the matrix element due to 
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Bz is r^Qji'^ = —jS'^^Bz, while away from the anti-crossing Therefore the second order contribution to the matrix 



orb 



-liBzhlX^- 



(35) 



Contrary to the case of electrically induced spin reso- 
nance, the oscillating magnetic field can induce transi- 
tions also at zero static magnetic field, as seen in Fig. [3};. 
However, the matrix element of the magnetic field is over- 
laid if the electric field in the x direction is present, since 
such electric field is much more efficient for the orbital 



resonance. 



(36) 



because it couples unperturbed states directly. 

If the electric field is oriented along y, it is much less 
effective, because the linear spin-orbit terms do not con- 
tribute in the first order. Here, for a non-zero matrix 
element in Eq.((34l). the perturbation Hi has to contain 
a term which is spin diagonal with spatial symmetry xy. 
The only such in Hi is the term originating in the first 
term of Hj^^, Eq. (fT7|) . After the rotation of the coordi- 
nate system this term is —{2'^c/Ibr) cos{2S)kxky, leading 
to the matrix element 

_ c- 7c „.,..w 4(g?"-gl") 
^ Vb - e^y^ 1 cos zoji r _ ^^0)2 _ (^lo _ ^oo)2 ■ 



(37) 



In small magnetic fields (< 1 T) this contribution domi- 
nates the matrix element compared to the contributions 

(2) 

from other parts in Hi, such as H^ , contributing in the 
second order. Note that there is no term with appropri- 
ate symmetry (spin diagonal, spatially xy) of Hi coming 
from a mixture of Hjj and Hjj^, making ^^^.^^ a specific 
effect due to the mixed cubic Dresselhaus and Bychkov- 
Rashba interactions. This example demonstrates the use- 
fulness of information about the symmetry contained in 
Eqs. (|15p - (|17p . By simple inspection of the symmetry we 
can tell immediately which term needs to be considered 
for a specific situation. 

The dependence of the matrix elements on the static 
magnetic field orientation 7 is shown in Fig. 3] c and d. 
The magnetic field matrix element is proportional to hf, 
see Eq. ([55]) . The direct couphng through the electric 
field along x is independent on 7. The matrix element 
of the electric field along y, as given in Eq. (|37p . is inde- 
pendent on 7 and can not be put to zero by changing the 
magnetic field orientation - as seen in Fig. 3):. However, 
there is some dependence to be seen and the dependence 
is striking for a different dots' orientation. The reason is 
that Eq. (j37p is the dominant contribution to the matrix 
element only up to a certain value of the static magnetic 

field - in higher fields the second order contribution from 

(2) 

H^ will dominate. Since there is already a visible de- 
pendence in Fig. 21;, we can estimate the crossover mag- 
netic field to be 1 Tesla, for our parameters. In Fig. HJi, 
the contribution of Eq. (|571) is zero exactly, since 6 = 45°. 



(2) 

element coming from H^ ' is seen. The possible depen- 
dence of the matrix element on 7 can decide whether the 
matrix element is induced by linear spin-orbit terms (de- 
pends on 7), or the mixed cubic- linear terms (does not 
depend on 7). This could be used as a detection for the 
presence of the cubic Dresselhaus term. Unless the elec- 
tric field is positioned exactly along y axis, no oscillating 
magnetic field infiuence or anisotropy can be observed 
due to high effectiveness of the electric field along x. 

After having analyzed means of control over the field 
matrix clement, or, in another words, Rabi frequency, we 
will now study the steady state solution of the density 
matrix. We will show that the Rabi frequency and de- 
coherence, which have been obtained in Refs. 'w.TT' from 
the observation of the decaying Rabi oscillations, can be 
obtained alternatively from the steady state current mea- 
surement. 



VI. RESONANT FIELD INFLUENCE IN THE 
STEADY STATE 

In this Section we are interested in the steady state 
solution of the density matrix, denoted by p and defined 
as the solution with constant occupations 



{dt + df)p.,,^0, Vt, 



(38) 



where the two contribution to the time derivative are 
those in Eqs. (O and ([9]). Even though it is not currently 
measurable in a single electron system, we include in our 
list of interesting steady state parameters the absorption. 



(39) 



defined as the energy gain of the electron due to the os- 
cillating field. 

After the decay of the Rabi oscillations, the system is 
in the steady state, where the occupations are constant. 
In this case the time derivative of the density matrix due 
to the oscillating field Eqs. ([9|) can be simplified to (see 
Ref. m for the derivation) 



df Paa = ~df pbb = 2{pbb - Paa)J, 



where the induced rate 

J = 



Iba 



4 A^+^l- 



(40) 



(41) 



The zero time derivative of the occupations in the 
steady state can be interpreted as a balance between two 
competing processes - relaxation [Eqs.©] which drives 
the system towards the thermodynamical equilibrium 
{pbb/Paa — ^ab/^ba) and osciUating field induced tran- 
sition [Eq. (HD)l ] equilibrating occupations of the resonant 
states (pbf, = Paa)- The effectiveness of the oscillating 
field in driving the system out of thermal equilibrium is 
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characterized by the induced rate J, Eq. pT|) . Going 
away from resonance the oscillating field is less effective 
in influencing the system, reflected by the (Lorentzian 
shape) decay of the induced rate. 

Our numerical strategy to obtain the steady state den- 
sity matrix p is as follows: We diagonalize the coupled 
dots electron Hamiltonian, Eq. ([T|)^ and compute the 
relaxation rates using Fermi's Golden rulci^ We choose 
a pair of resonant states, {a, b}, and after evaluating 
^Iba we find the induced rate according to Eq. ([1T|) . Fi- 
nally, we find the steady state density matrix by solving 
the set of linear equations defined by Eq. (|38p . A differ- 
ent method, with the oscillating field treated exactly, was 
used for single dot in intense oscillating fields;^ three or- 
ders of magnitude larger than the fields considered here. 

We can analytically reproduce the numerical results by 
the two state approximation discussed in the above. The 
physics is then characterized by the number 



"^0 — ^ba'-^\'^='-Jba — 1^ 



1 



6a I 



476a EbQ 



(42) 



which is the induced rate at the resonance, measured 
in the units of the relaxation rate between the resonant 
states. 

Two limits can be identified, according to Jq. If the 
induced rate dominates the relaxation, Jq 1, the oc- 
cupations of the two resonant states are close to being 
equal, while if Jq <C 1, the system is close to the ther- 
mal equilibrium. The interpretation of 2 J as the electron 
outscattering rate due to the oscillating field, as it follows 
from Eqs. (HO]), is reassured by the result form the ab- 
sorption. We expect the absorption to be proportional 
to a transition rate from the excited state to the ground 
state times the energy dissipated at this transition. If 
Jg <C 1 the transition rate is 2J. In the opposite limit, 
Jg 1, the outscattering due to the oscillating field is 
strong and the transition rate for the dissipation is lim- 
ited by the relaxation rate. The frequency full widths at 
half maximum (FWHM) also differ for the two limits ~ 
see Tab. HI] for analytical results. 

Figure [S] presents our numerical results for induced 
rate, excited population width, and decoherence as func- 
tions of the tunneling energy for the spin and orbital 
resonance. Both resonances are in the regime of Jg ^ I, 
where the decoherence is revealed by the FWHM of the 
induced rate, see Tab. [TTl while the relaxation rate can 
be obtained if both the induced rate at resonance and 
FWHM of the excited population are known, too. Due 
to Eq. (fTU)) . the relaxation rate is indiscernible from the 
decoherence in the figure and Jg can be directly deter- 
mined. For the spin resonance Jq varies between 10^ 
and 10^^ - the limit expressions in Tab. [Hlare then ex- 
act with this precision. The upward dips in FWHM and 
the decoherence rate are due to the anti-crossing of the 
spin and orbital stateSf^^ It is interesting that the in- 
duced rate is not influenced by the anti-crossing. This is 
because both the square of the matrix element and the 
decoherence (equal to the relaxation) in Eq. (PT|) depend 



a 



Pbb 
J 
W 



Pbb 
Pbb 
W 

w 



steady state 



2J+r6„(l + r) 



EbaJ 



l + r + 2.J/r6„ 



at resonance 



2JJ + 1 + T 



\ + T + 2.J, 



limit 



JL » 1 
Jo « 1 
Jo » 1 
Jo" « 1 



FWHM {Sul 



12} 



8J6(l + J6)+4^(l+^+3Jo) 3 



-t(1+t+3JJ) 

Ha 

4(l+r+2JJ2 2 
1 + T- Iba 



Iba 



at resonance 



1/2- (l-r)/2J^ 
^,+JE{l-r)/{l + T)- 

Eba^bai'^ - r) 

2gft„J'-(l-r)/(l + r) 



FWHM {Suol 



111 



2\Q.^aflbalVba{\-'iT) 

'A^baflbalVbaiX + r) 
476a 



TABLE H: (a) Steady state, value at resonance, and frequency 
full width at half maximum (FWHM) 50)1/2 squared for the 
excited state occupation p;,;,, the induced rate J, and absorp- 
tion W . Note that the FWHM of the excited population 
is defined only if the temperature is low enough such that 
Jo ^ ''"(1 + ''")/(l — 3r). (b) The value at the resonance, and 
frequency full width at half maximum of the excited popula- 
tion and absorption in the two limits. 



on the anti-crossing in the same way and the contribu- 
tions cancel. Also note that the rates characterizing the 
oscillating field are very different in the transient and the 
steady state regime. While the steady state characteris- 
tic rate J is ~ 10^'^ s~^, looking at Fig. [SJd one can see 
that the Rabi frequency for the same parameters is only 
^ 10« s-i. 

Compared to the spin resonance, the orbital resonance 
is much less sensitive to the anti-crossing, since only in 
a very narrow region at the anti-crossing the relaxation 
rate acquires a factor of one half.^^ One also sees that Jq 
is smaller, meaning we are closer to the regime of Jg < 1 
which can be reached by lowering the amplitude of the 
oscillating electric field. In that regime, the decoherence 
can be obtained from the FWHM of the excited popula- 
tion or from the induced rate. 

We finish this section by summarizing that after iden- 
tifying the appropriate regime of high or low induced rate 
one can obtain the decoherence and Rabi frequency us- 
ing expressions from TabHIl provided one can measure 
the induced rate J and the excited state population pj,^ 
(and their full widths). In turn, these two can be mea- 
sured if the dot is connected to leads and the current 
flows through the dot, as shown theoretically in Ref. [sil . 
In Ref. i23| it is shown even on a simpler model that the 
measurement can be done by changing the coupling be- 
tween the dot and the leads. Namely, for small coupling 
the current is proportional to the excited state popula- 
tion, while for large coupling the current measures the 
induced rate. 



VII. CONCLUSIONS 

We have studied electrically and magnetically induced 
spin and orbital resonance of a single electron confined 
in coupled lateral quantum dots. We have taken into ac- 
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FIG. 5: Calculated induced rate J at resonance (solid), de- 
coherence 7i,a (dashed), and the FWHM of the excited popu- 
lation (dot-dashed) as functions of the ratio of the tunneling 
energy 5Et and the confinement energy Eq for (a) spin reso- 
nance and (b) orbital resonance. The static in-plane magnetic 
field is S = 1 T. If the solid line is above (under) the dashed 
one, it means that JJ > 1 (JJ < 1). The dots are oriented 
along [100], while the static magnetic field lies along [010]. 



count the relaxation and decoherence due to an acoustic 
phonon environment, with the rates computed by Fermi's 
Golden rule. Resonant oscillating electromagnetic fields 
are capable to induce transitions between electron eigen- 
states. We have focused on the oscillating field matrix 
elements, equal to the Rabi frequency, for spin and or- 
bital resonance. 

We have given an effective spin-orbit Hamiltonian 
which allows to quantify the spin-orbit influence on the 
matrix element using symmetry considerations. Specif- 
ically, for electrically induced spin resonance, we have 
shown how the spin-orbit anisotropy allows to control 
the matrix element by the strength and orientation of 



the static magnetic field. These conclusions give hints for 
optimal quantum dot configurations for the case when: 
(i) the spin is manipulated by an oscillating electric field, 
whereas its influence is desired to be maximized and (ii) 
the spin manipulated by an oscillating magnetic fleld, 
when the effect of the electric field on the spin is desired 
to be minimized. Connecting with our previous work, we 
have found that the easy passage provides not only long 
spin relaxation time, but also stability against electric 
field disturbances, making it a suitable arrangement for 
spin qubit realization. 

In a double dot, the electric field is most effective in 
spin manipulation if it lies along the dots' axis and the 
matrix element is strongly influenced by the anti-crossing 
(a crossing of different spin states lifted by spin-orbit in- 
teractions). An important feature is that the electric 
field is less effective if the magnitude of the static mag- 
netic field is lowered. Oscillating electric field of order 
of 1000 V/m can easily be more effective than oscillating 
magnetic field of 1 mT if the static magnetic field is of 
order of Tesla. For these parameters in a GaAs quantum 
dot the Rabi frequency of 1 GHz is achievable for the 
spin manipulation using an electric field. 

In the last part we studied the influence of the resonant 
fields in the steady state. We proposed the induced rate 
as a single characteristic paraiiic;t(^r. We have analyzed 
steady state occupations, induced rate, and absorption 
and their full widths for both spin and orbital resonance 
and used those results to show how to obtain decoherence 
and Rabi frequency from these steady state characteris- 
tics. In turn, these characteristics can be obtained from 
a steady state current measurement. 
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